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ABSTRACT 


The objective of this study is to introduce three estimators of parameter Weibull distribution, depending on 
standard Bayes estimation and modified Bayes estimation. Then, using a simulation study (MATLAB program), to find the 
best one base on MPE and MSE. 
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INTRODUCTION 

The theorem of Bayes provides a solution on how to learn from data. Related to survival function and by using 
Bayes estimator, Elli and Rao (1986), estimated the shape and scale parameters of the Weibull distribution, by assuming a 
weighted squared error loss function. They minimized the corresponding expected loss, with respect to a given posterior 
distribution. Sinha & Sloan (1988) obtained, Bayes estimator of three parameters Weibull distribution and compared the 
posterior standard deviation estimates, with the corresponding asymptotic standard deviation of their maximum likelihood 
counterparts and numerical examples are given. In 2002, Klaus Felsenstein developed Bayesian procedures, for vague 
data. These data were assumed to be vague in the sense that the likelihood is a mixture of the model distribution with error 
distribution. An extension of Jeffery prior information with square error, loss function in exponential distribution was 
studied by Al-Kutubi (2005). In this paper, Alkutubi (2009) proposed an extension of Jeffery prior information, with a new 
loss function and then compare it with standard Bayes, to find the best. In this paper, we will introduce three estimators to 
the parameters Weibull distribution, using standard Bayes estimation and extension loss function. We will then compare 
between them by a simulation study, to find the best one base on MPE and MSE. 

MATERIALS AND METHODS 

Let tj, t 2 ,..., t n be the life time of a random sample of size fl with distribution function and probability density 
function. In the Weibull case, we assumed that the probability density function of the lifetime is given by 



To obtain Bayes estimator, the following steps are needed [2], [3]. 


A number of fl items put to the test and the lifetimes of this random samples are recorded with the probability 
density function f(t, 6, c).The life time ability density function f(t, 0, c) are regarded as a conditional probability density 

function f(t\9, c ), where the marginal probability density function of 6 is given by §(y) , the Jeffery prior information. 
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We find Jeffery prior by taking ,, (0 ) c* 1 (0 ) , where fisher information/(0) is given by 


, w = _„ E (£i^£M2) 



Then 1(9) — so 


the Jeffery prior information is given byg(9) — k^J I (9) — k 



The joint probability density function is: 


e n+1 ve7 


exp 



The marginal probability density function of (t L ,... t n , 6, c) is given by: = 


fccVn Zn-lN | 

(sto n \—)' 


Then posterior distribution is, 


cffiti )’ 1 



e n+1 (^) ! 


By using squared error loss function= 


, we can obtain the Risk function, such that 


R(9,9) =EL(9,9) = f*c(9-9) 2 n(9\t 1 ,...t n )d9 

, dR(e,e) _ , „ . ■ -7r Y.t t 

Let—= 0, then Bayes estimator is 9 1 - ^ n _ 1 

By using new extension loss function [1], we can get new estimator of parameter Weibull distribution, such that 

R(0’9) = %9‘(9-9) 2 7rC@ I t-L,... t n )d9 Solving equation in above to get the second Bayes estimator 

a — £L 
° 2 ~ (?1 i 

By c=2, The another new estimator of parameter Weibull distribution, such that 

R(8’9)=£’9 2 (9-9) 2 7t(0 I ... t n )d9 Solving the above equation to get the second Bayes estimator 



RESULTS 


In a simulation study, we have chosen n=50, 75, 100 to represent small, moderate and large sample size, 
several values of parameter #=0.5, 1, 1.5, and values of the new loss function c = 0.4, 0.8. The number of replications 
used was R=1000. The simulation program was written by using the Matlab program. After the parameter was estimated, 
mean square error (MSE) and mean percentage error (MPE) were calculated, to compare the methods of estimation, where 
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The results of the simulation study are summarized and tabulated in Table 1 and Table 2 for the MSE and the 
MPE of the three estimators for all sample sizes and lvalues respectively. It is obvious from these tables, 
a Bayes estimator with the Jeffery prior information, 6^is the best estimator. But 0 3 in most cases have the largest MSE 
and MPE 

Table 1: The Ordering of the Estimators with Respect to MSE 


Size 

e 


Theta Flat 



0, 

0 3 


0.5 

0.4 

0.100 

0.135 

0.211 


0.8 

0.121 

0.133 

0.211 

50 

1 

0.4 

0.182 

0.209 

0.229 

0.8 

0.189 

0.209 

0.229 


1.5 

0.4 

0.214 

0.231 

0.233 


0.8 

0.223 

0.235 

0.233 


0.5 

0.4 

0.125 

0.131 

0.139 


0.8 

0.125 

0.132 

0.132 

75 

1 

0.4 

0.121 

0.122 

0.122 

0.8 

0.123 

0.123 

0.123 


1.5 

0.4 

0.122 

0.125 

0.130 


0.8 

0.126 

0.126 

0.133 


0.5 

0.4 

0.091 

0.093 

0.094 


0.8 

0.091 

0.093 

0.094 

100 

1 

0.4 

0.088 

0.089 

0.089 

0.8 

0.088 

0.089 

0.089 


1.5 

0.4 

0.088 

0.088 

0.089 


0.8 

0.088 

0.088 

0.090 


Table 2: The Ordering of the Estimators with Respect to MPE 


Size 

0 


Theta Hat 


0, 

07 

07 


0.5 

0.4 

0.010 

0.013 

0.021 


0.8 

0.010 

0.013 

0.021 

50 

1 

0.4 

0.012 

0.020 

0.022 

0.8 

0.012 

0.020 

0.022 


1.5 

0.4 

0.013 

0.013 

0.024 


0.8 

0.013 

0.013 

0.024 


0.5 

0.4 

0.012 

0.012 

0.013 


0.8 

0.012 

0.012 

0.014 


1 

0.4 

0.011 

0.013 

0.014 

J 

0.8 

0.011 

0.013 

0.014 


1.5 

0.4 

0.010 

0.013 

0.013 


0.8 

0.010 

0.013 

0.012 


0.5 

0.4 

0.006 

0.007 

0.009 


0.8 

0.006 

0.007 

0.009 

100 

1 

0.4 

0.005 

0.006 

0.007 

0.8 

0.004 

0.004 

0.007 


1.5 

0.4 

0.002 

0.003 

0.005 


0.8 

0.002 

0.003 

0.005 


DISCUSSION 


The results of the simulation study are summarized and tabulated in Table 1 and Table 2 for the MSE and the 
MPE of the three estimators for all sample sizes, 0and C values respectively. The order of the estimator is from the best 
(smallest MSE) to the worst (largest MSE). It is obvious from these tables, a Bayes estimator with Jeffery prior information 
and standard loss function, 6 1 is the best estimator. In most of the cases, it is apparent that, the Bayes estimator with Jeffery 
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prior information and extension new loss function d 2 i s the next best estimator. But, standard Bayes estimator with new 
loss function § 3 in most cases has the largest MSE and MPE. 

CONCLUSIONS 

The Bayes estimator, with Jeffery prior information 0 1 is the best estimator, when compared it with the other 
estimators. We can easily conclude that, MSE and MPE of Bayes estimators decrease with an increase of sample size. 
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